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ABSTRACT 

> 

^ ■ We compare the results of Eulerian hydrodynamic simulations of cluster formation 

against virial scaling relations between four bulk quantities: the cluster mass, the dark 
matter velocity dispersion, the gas temperature and the cluster luminosity. The comparison 
is made for a large number of clusters at a range of redshifts in three different cosmological 
models (CHDM, CDM and OCDM). We find that the analytic formulae provide a good 
description of the relations between three of the four numerical quantities. The fourth 
(luminosity) also agrees once we introduce a procedure to correct for the fixed numerical 
resolution. We also compute the normalizations for the virial relations and compare exten- 
sively to the existing literature, finding remarkably good agreement. The Press-Schechter 
prescription is calibrated with the simulations, again finding results consistent with other 
authors. We also examine related issues such as the size of the scatter in the virial relations, 
the effect of metallicity with a fixed pass-band, and the structure of the halos. All of this is 
^ ■ done in order to establish a firm groundwork for the use of clusters as cosmological probes. 

^ . Implications for the models are briefly discussed. 

Subject headings: galaxies: clusters, X-rays: galaxies, methods: numerical 
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1. Introduction 

The statistics of X-ray clusters can serve as an excellent 
probe of cosmology. The luminosity function of clusters is 
defined in a straightforward manner both observationally 
( [Henry et al. 1992fc |Ebcling ct al. 19951; |Burns et al. 1996| ) 



and numerically ( Cen 1992| ; Bryan ct al. 19941: ) , although 
as we shall see this is a computationally difficult task. 



The temperature function is similarly obtainable (Henry 



& Arnaud 1991 ; David et al. 1993 ), however, here it is the 
observational data which are more challenging to obtain. 

One uncertainty in performing this comparison is the 
possibility of systematic errors in the theoretically de- 
rived cluster temperatures and luminosities. This has 
been investigated in a number of ways, primarily by test- 
ing individual methods developed against simplified prob- 
lems with known solutions. So, the Lagrangian smoothed 
particle hydrodynamics (SPH) method combined with a 
P 3 M code for computing the gravitational interaction has 
been examined by Evrard (1988), Couchman, Thomas 
& Pearce (1995) and others; while SPH with a gravita- 
tional tree code has been tested in a number of papers 



Another way to check the results of simulations (and 
the route taken in this paper) is to compare them against 
the predictions of approximate analytic models. Al- 
though agreement does not guarantee correctness (as 
both methods are only approximations to the full solu- 
tion), concordance would increase our confidence in both 
methods. Also, simple analytic arguments may only iden- 
tify a scaling property between quantities without speci- 
fying a normalization, which can however, be fixed by nu- 
merical simulation (or by further assumptions). This de- 
scribes the scaling laws that come from considering clus- 
ters as spherical clouds of gas in hydrostatic equilibrium. 
Navarro, Frenk & White (1995) (hereafter NFW) recently 
compared the results of six clusters simulated with SPH 
in a Cold Dark Matter (CDM) scenario against these scal- 
ing relations (at z = 0). They find good agreement over 
a wide range of luminosity, mass and temperature, but 
claim that clusters from Eulcrian simulations (such as in 



Kang et al. 1994a and Bryan et al. 1994a) do not 



(Hcrnquist & Katz 1989; Katz, Weinberg & Hcrnquist 



Another analytic method is that initially described by 
Press & Schechter (1974) which predicts the mass distri- 
bution of collapsed objects. There have been a number 
of comparisons between its predictions and the results of 



L9 96| |Navarro fe White 1993]). A novel modification of N _ body simulations flEfstathiou et jOggf ; |Bond et al 



SPH was described by Shapiro et al. (lyyo), along with a 
number of comparisons against known results. For Eule- 



1991 ; Lacey fc Cole 1996 ). Using the scaling results, this 



ria* codes, method papers with some tests include: Cen 
(1992) for a first-order grid based code; Ryu et al. (1993) 
for a total variation diminishing method; Anninos, Nor- 
man & Clarke (1994) for a two-level nested grid scheme, 
and Bryan et al. (1995) for a code based on the piecewise 
parabolic method. 

However, a number of recent papers have taken up the 
question of accuracy and consistency with a cosmological 
context as their primary focus. For Eulerian codes (with 
fixed comoving spatial resolution) , this issue has been ad- 
dressed by Anninos & Norman (1996) (hereafter AN96) 
who performed simulations with a two level hierarchical 
method and simulated the formation of a moderately rich 
cluster with five different resolutions ranging from 1600 
to 100 kpc per cell. They found that although the tem- 
perature, mass and velocity dispersion of a cluster were 
reasonably well determined with the lower-resolution sim- 
ulations, the total luminosity had not converged even for 
the highest resolution run. The bolometric luminosity 
(for this single cluster) behaved as 

where Ar is the spatial resolution of the simulation and 
a = 1.17. 



theory can be extended to produce the temperature (Eke 



Cole & Frenk 1996) and luminosity distribution functions. 



Since this is one of the most widely used constraints on the 
amplitude of mass fluctuations, it is important to check 
its validity. 

In this paper, we make a detailed comparison between 
simulation results and the adiabatic scaling laws as well 
as the Press-Schechter formalism with extensions. This 
allows us to gauge the accuracy and consistency of both 
methods, leading to firmer conclusions regarding the vi- 
ability of the cosmology modelled. The paper is laid out 
as follows. In section |^, we review the scaling relations, 
including a modification to take into account the finite 
resolution of Eulerian codes. We then compare these 
to the results of CDM and Cold plus Hot Dark Matter 
(CHDM) simulations at a variety of redshifts. In sec- 
tion U, we examine the mass, temperature and luminosity 
distribution functions, including the effects of finite band- 
pass and line emission. These are compared against the 
Press-Schechter plus scaling theory (extended to include 
the additional complications in the luminosity function). 
In section || we briefly examine the profiles of temper- 
ature and velocity dispersion to assess the accuracy of 
the isothermal models assumed in extending the Press- 
Schechter work. Finally, in section |^, we discuss our re- 
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suits and comment on the viability of the models simu- 
lated. 



2. Scaling Relations 



The value of A c is taken from the solution to the 
collapse of a spherical top-hot perturbation under the 
assumption that the cluster has just virialized ( |Pccblcs 
1980). Its value is 18-7T 2 for a critical universe but has a 



Here we review the scaling relations between cluster 
bulk properties through the assumption of a pressure sup- 



poncea isocncrmai sphere 

and oiie-dimeiisioiial collisioiiless velocity dispersion a of 

the dark matter particles. The assumption of a specific 
density profile (here the isothermal sphere) is not required 
to obtain the scaling behaviour, but is needed to deter- 
mine the constant of proportionality between the given 
quantities. 

These relations were used by Kaiser (1986) to describe 
the evolution of 'characteristic' quantities, largely driven 
by the non-linear mass (M„;), defined via equation ( |l9| ) 
below, as well as to derive relations between distribution 
functions at different epochs. We do not explicitly test 
these because they are uniquely specified by the non- 
linear mass (which we do examine) and the scaling re- 
lations discussed below. 

2.1. Scaling review and normalization 

In the isothermal distribution function, the density is 



dependence on cosmology through the parameter fi(z) = 
ilo(l + z) 3 j E{z) 2 . We have calculated this for the cases 
(|Lacey fc Cole 1993| ) and n R = ( [Eke, Cole fc 



properties mrougn me assumption oi a pressure sup- Q A = Q QLacey fc Cole 1993| ) and ft 
ijted isothermal sphere for both the gas temperature T Frcnk 19961) fitting the results with - 
d oiie-dimeiisioiial collisioiiless velocity dispersion a of 



relate d to the velocity dispersion ( Binney 
1987b : 



Trcmainc 



p(r) 



(2) 



2-nGr 2 ' 

If we define r V i r as the radius of a spherical volume within 
which the mean density is A c times the critical density 
at that redshift (M = 4:irr 3 ir p crit A c /3), then there is a 
relation between the virial mass and the one-dimensional 
velocity dispersion: 



cr = M 1 / 3 (H 2 {z)A c G 2 /l6) 



l/G 



(3) 
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10 15 M< 



o 



1/3 



7i 2 A c £ 2 ) 1/6 km/s (4) 



In the second line we have introduced a factor f a which 
will be used to match the normalization from the simu- 
lations. The redshift-dependent Hubble constant can be 
written as H(z) — 100 h E(z) km/s/Mpc with the func- 
tion E 2 = Sl (l + z) 3 + Ojcj (1 + z) 2 + fl A dependent on 
three contributions: 



n 



8irGpo 



9.r 



1 



(H {) Rf 



A 
3H 2 ' 



(5) 



Here, po is the non-relativistic matter density, R is the 
radius of curvature and A is the cosmological constant. 



A c = 18tt 2 + 82x - 39x 2 for fl R = (6) 
A c = 18tt 2 + 60x - 32x 2 for il A = 

where x = fl(z) — 1. These are accurate to 1% in the 
range tt(z) = 0.1-1. 

If the distribution of the baryonic gas is also isothermal 
wc can define a ratio of the 'temperature' of the collision- 
less material (T CT — pm p a 2 /k) to the gas temperature: 







[xrripCT 
kT 



(7) 



We take /i = 0.59. Given equations (j^) and (fjj), the 
relation between temperature and mass is: 



kT = 



GM 2 / 3 pm p { H 2 (z)A 
1.39/t 



AI 



2G 

2/3 



1/3 



(8) 



10 15 M m 



(/i 2 A c £ 2 ) 1/3 keV (9) 



Since we will test the relations ||, and || against the 
numerical simulations separately, we also add a normal- 
ization factor, fx, and set j3 = 1 for this last equation. 

We can easily find the scaling behaviour of a cluster's 
X-ray luminosity by assuming bolometric Bremsstrahlung 
emission and ignoring the temperature dependence of the 
Gaunt factor (e.g. |Spitzer 197g| ): L boi cx MpT 1 / 2 . 

We could compute the luminosity by using the isother- 
mal sphere approximation, however this is either infinite, 
if there is no rollover in density as r — * 0, or we must ar- 
bitrarily select a core radius. Instead, we will assume the 
normalization found by NFW (adjusted slightly to match 
our definition of the virial mass): 



L. 



hot 



1.3x10 



45 



M 



4/3 



(h 2 A r E 2 



7/6 ( rib 



n) erg/s ' 

(10) 

To obtain this we have used the redshift dependence of 
the critical density, the temperature-mass relation and 
multiplied by the mean baryon fraction. 
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Other scaling laws can easily be derived from these; we 
will need one more: 

i — »i^)Wr(|) 2 « g / 8 . 

(id 

We note that cosmology enters these relations only with 
the combination of parameters h 2 A c E(z) 2 , which comes 
from the relation between the cluster's mass and mean 
density. The redshift variation comes mostly from E(z) 
which is equal to (1 + z) 3 / 2 for an Einstein-de Sitter uni- 
verse. 

2.2. Resolution effects 

While we do not expect the first three bulk proper- 
ties (mass, temperature and velocity dispersion) to be 
strongly dependent on the numerical resolution, the X- 
ray emissivity is very sensitive to the density profile of 
the cluster since it originates primarily from the central 



region ot the cluster, in order to examine trie expected 
luminosity behaviour with resolution, we take the density 
profile fit from high-resolution N-body simulations found 
in Navarro, Frenk & White (1996): 



So/c 3 



P( r ) _ 

Peril x(l/c + x) 



2 ' 



(12) 



where x = r/r V i r . For clusters, a reasonable parameter 
choice is c = 5, although our results will not be sensitive 
to small changes in this value. This is a fit to the colli- 
sionless component and not the baryonic gas (as would be 
more appropriate), however the difference only becomes 
important at r/r V i r < 0.04 (c.f. Navarro, Frenk & White 



1995; Bryan & Norman 1997a) which will not substan- 



tially affect our results. 

In order to approximate the effect of finite numerical 
resolution we filter the density distribution: 

P'(r,r sm ) = — „2 J" '- . (13) 



The smoothing kernel is a Gaussian (Wg(i) = e~ x2 l 2 ) 
and r sm is the smoothing radius. Using symmetry, this 
integral can be partially computed analytically: 



1 



P = 



27TTT 



M J 



-(r'-r) 2 -jr' + r) 2 



p{r )r e 



(14) 

In Figure [j]a, we show the effects of various smoothing 
radii to the adopted density profile. The mass is con- 
served by the smoothing process in the sense that the 



total mass within some fixed fraction r/r sm ^ 1 is con- 
stant as r sm changes (the total mass of equation ( |l2| ) 
does not converge as r — > oo). This results in mass being 
transferred from small to large radii, causing the profile 
to steepen outside the smoothing radius, another feature 
observed by AN96. Assuming the temperature profile 
does not change, the bolometric free-free luminosity can 
be computed as a function of r sm /r V i r and the result is 
shown in Figure |l|b. Also shown is the scaling result found 
by AN96 and given by equation ([!]) , which they found to 
be approximately valid in the range 0.1 > r sm jr vir > 1. 
The agreement is good enough that we will adopt their 
value (a = 1.17), although it should be kept in mind that 
this power-law behaviour is only approximately correct. 

If all density profiles can be scaled by the virial radius 
and critical density to agree with equation ([l2]) regardless 
of mass (at least in a statistical sense), then the smoothed 
luminosity will scale in general as Ltoi oc (r V j r /r srn ) a . R e- 



cent work ( [Navarro, Frenk fc White 1996| ; |Colc fc Lacey 



1996 ) indicates that, in fact, this is not true and lower 
mass objects form earlier and have steeper profiles. Nev- 
ertheless it is a good approximation over the range of 
masses examined here and we will adopt it. 

Using the relation between the virial radius and virial 
mass, we can include the effects of numerical resolution 
in equation (|l0|): 



hoi 



5.8 x 10 



45 



M 



10 15 M Q 

/ 100 kpc 
I Ax 



l + o 

3 



(h 2 A c E(z) 2 ) 



(1 + z) a erg/s, 



(15) 



Here we have written the resolution in terms of the co- 
moving cell size Ax and normalized by assuming that 
Ax = 100 kpc will correctly reproduce the observed core 
radius of a 10 15 M Q cluster (c.f. AN96). The resulting 
L' bol -T relation is, 



a 



hoi 



2.5 x 10 



15 



kT/f T \ ■ 
1.0 keV/ 



(h 2 A c E(z) 2 



n b \ 2 /ioo k P c\ Q 



n 



J (1 + z) a erg/s 



(16) 



2.3. Comparison to simulations 

Our numerical techniques have been described in detail 



elsewhere (Bryan et al. 1995), so we only briefly summer- 
ize them here. The dark component was modelled by the 
particle-mesh method with the Poisson equation solved 
through an FFT. The adiabatic equations of gas dynam- 
ics are solved with a modified version of the piecewise 
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Fig. 1. — (a) The density profile in equation ( |l2| ) 
smoothed by a Gaussian filter with three different val- 
ues of the smoothing radius, (b) The relative bolomct- 
ric bremsstrahlung luminosity of the smoothed density 
profile as a function of smoothing radius. The dashed 
line shows the result found by direct numerical simula- 
tion (AN96). 



parabolic method. Both techniques have a resolution of 
a few cell widths. Radiative cooling is not included as we 
have insufficient resolution to properly follow the cooling 
structures; however, for large clusters, which have very 
long cooling times over most of their volume, this may be 
a reasonable approximation. 

The four simulations we use in this paper are summa- 
rized in Table § The results from the CHDM256 model 
are sufficiently similar to CHDM512 that we mostly focus 
on the first three models. The resolution, as measured by 
cell width, is similar for the three 85ft. -1 Mpc boxes, while 
substantially better for the CHDM512 model, giving us 
some leverage in studying the effects of resolution. All 
have power spectrum normalization (indicated here by 
the as, the size of mass fluctuations in spheres of radius 
8ft" 1 Mpc) which are approximately in agreement the the 
amplitude of fluctuations in the cosmic background radi- 
ation on large scales. Some of these simulations and their 



initializations have been described elsewhere ( Bryan et al 
19944 |Bryan et al. 1994b| ). 



To identify clusters we adopt the spherical overden- 
sity algorithm described in Lacey & Cole (1996). Peaks 
in the density distribution are identified as cluster cen- 
ters; spheres are grown around each point until the mean 
interior density reaches A c . The center-of-mass inside 
the region is computed and this points is used to grow 
a new sphere. The procedure is iterated to convergence 
and then all clusters are checked for overlap (defined as 
a center-to-center distance less than three-quarters of the 
sum of their virial radii), with the less massive cluster be- 
ing removed. The scaling relations are robust to changes 
in what we do with overlapping clusters. 

In the following sections we examine each scaling law 
in detail and compare to other work where possible. 

2.3.1. The M-a relation 

In Figure |^, the mass-weighted one-dimensional veloc- 
ity dispersions and virial masses are plotted for clusters 
identified in our three primary models, at three different 
redshifts: z = 0, 0.5 and 1 (except for the CHDM512 
simulation, for which the z = 0.5 output was corrupted 
so z = 0.3 was used instead). The velocity dispersion in- 
cludes both cold and hot collisionless components, where 
appropriate, and is computed after removing the center- 
of-mass velocity. There are substantially more clusters in 
the CDM270 simulation because of its higher normaliza- 
tion (the CHDM512 run also used a somewhat smaller 
volume) . 

We also show the predicted virial relation between 
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Fig. 2. — The virial mass and one-dimensional velocity dispersion of clusters identified in the three models for three 
different redshifts: z = 0, 0.5 (0.3 for CHDM512) and 1 along with the scaling relation from equation (|j) with f a given 
by Table |. 



(i 



Table 1 
Simulation parameters 



designation 


^cold 


n hot 


^baryon 


h 


m v (eV) 


C8 


^cell 


Npart 


L box (ft-lMpc) 


CDM270 


0.94 


0.0 


0.06 


0.5 





1.05 


270 3 


135 3 


85 


CHDM512 


0.725 


0.2 


0.075 


0.5 


2 x 2.3 


0.7 


512 3 


3 x 256 3 


50 


OCDM270 


0.34 


0.0 


0.06 


0.65 





0.75 


256 3 


128 3 


85 


CHDM256 


0.6 


0.3 


0.1 


0.5 


7.0 


0.67 


256 3 


3 x 128 3 


85 



these two quantities, given by equation ^ . The normal- 
ization, specified by / CT , is determined separately for each 
model by fitting the 25 most massive clusters at z = 0; 
the numerical values are listed in Table ||, along with their 
standard deviations. The results are almost completely 
insensitive to small (factor of two) changes in the number 
of clusters used. The same normalization is used at higher 
redshifts. For the most massive — and best resolved — 
clusters, the agreement in both slope and redshift vari- 
ation is quite good, with remarkably little scatter. For 
less well-resolved clusters, those with a below 600 km/s, 
the slope flattens. This is most likely due to resolution 
effects, as the virial radius drops to just a few cell widths. 

To examine any systematic variation with resolution, 
the normalization from each model is plotted against the 
cell size of that simulation in Figure |3| (here we include the 
CHDM256 model). Shown as an error bar is the scatter, 
parameterized by each sample's standard deviation. Al- 
though not very significant, the best fitting straight line 
shows a modest increase in f a with increasing resolution. 
This is in the expected direction since improved force res- 
olution requires higher velocity dispersions to maintain 
equilibrium. 

In Table we also list the normalizations determined 
by other authors, using different numerical techniques. 
Evrard (| vrard 1991] , E91) performed 22 simulations of 
isolated clusters in an Einstein-de Sitter universe with 
a cold dark matter spectrum and obtained a somewhat 
larger normalization. The six (£1 = 1) clusters of NFW, 
using a TREESPH code, also had a slightly larger nor- 
malization. They used a marginally different virial over- 
density, however when fitting their results to equation (^) , 
we adopt their value of A c = 200. The effect of increas- 
ing this parameter is to sample a slightly smaller volume, 
however since the velocity-dispersion profile is nearly flat 
in this range (c.f. section^), the effect is very small. The 
P 3 M simulations of Crone and Geller (1995, CG95) re- 
sulted in a mean f a — 0.91 (after converting from their 



1.00 



0.90 



0.80 



0.70 



0.60 



• velocity dispersion 
□ temperature 



100 200 300 400 500 600 700 800 
Ax (kpc) 

Fig. 3. — Normalizations for the M-T and M-a virial 
relations (as determined by numerical simulation) as a 
function of resolution, measured by cell size in kpc at 
z — 0. The f a points have been slightly offset in Ax for 
clarity. Lines are the best-fit linear relation. 
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Table 2 
vlrial fitting parameters 



source 


u 


Act 


fx 


AT 


^clust 


CDM270 


0.85 


0.05 


0.79 


0.10 


25 


CHDM512 


0.89 


0.03 


0.75 


0.14 


25 


OCDM256 


0.85 


0.04 


0.78 


0.05 


25 


CHDM256 


0.82 


0.04 


0.75 


0.05 


25 


E91 


0.98 


0.02 


0.81 


0.05 


22 


NFW 


0.96 


0.04 


0.86 


0.11 


6 


CG95 


0.91 








~50 


CL96 


~ 1 








~100 


EMN96 






0.92 


~0.05 


58 



notation) for well resolved clusters, with A c = 300. This 
also found that the smallest clusters deviated slightly 
from the virial power-law, ascribing the effect, as we do, 
to insufficient numerical resolution. Their scatter appears 
similar to that found here. Finally, the P 3 M simulations 
of Cole & Lacey (1996, CL96) produced similar results. 
They did not fit the virial relation, but their Figure 3 in- 
dicates both that the virial relation is a reasonable fit and 
that there is not substantial variation from a power law 
over the range of masses examined. Jing & Fang (1994) 
argued for f a = 0.8-0.9, but this was based on combined 
N-body and Press-Schechter comparisons, rather than di- 
rectly from the M — a relation. Finally, a number of other 
authors have studied this issue but often in ways which 
make it difficult to compare directly to the virial formal- 
ism. For example, Walter & Klypin (1996) used a fixed 
cluster radius rather than one based on the mean over- 



isothermal virial value: 
a 2 



2c 



ln(l 



dx 



1 



y 2 dy 



1 



1 + In ( 1 



(17) 



der sity. 



For the range of masses and cosmologies studied here, c 
varies from 5-10, and the resulting f a = 0.92-0.97. 

Thus we have strong evidence from a variety of sources 
that the virial relation describes clusters fairly well, al- 
though the normalization may need to be slightly modi- 
fied, with a value around f a ~ 0.9-0.95 in good agreement 
with a wide variety of numerical methods. Unfortunately, 
there is likely to be a systematic bias between the veloc- 
ity dispersion of the dark matter and that of the galaxies 
flCarlberg, Couchman fc Thomas 1990 ; Summers, Davis 
fc Evrard 1995| ) which we cannot address with these sim 



In section 2.2, we mentioned that high-resolution N- 
body simulations indicate that the dark matter density 
profile does not follow the isothermal r~ 2 profile but is 
instead given by equation ([l2|), which has two parame- 
ters. If one (So) is set by the requirement that the mean 
density inside the viri al radius is A c , the other is stro ngly 
correlated with mass ( Navarro, Frenk fc White 1996 ), so 
the mean radial profile is completely specified by the mass 
(or, more accurately, the redshift of formation). The im- 
plied value of f a (i.e. the correction as compared to the 
isothermal profile), can be computed using Jeans equa- 
tion (p~ 1 d[pa(r)' 2 }/dr = —d<f>/dr; we assume the velocity 
dispersion is isotropic for simplicity). It is the ratio of 
the mass weighted velocity dispersion (squared) to the 



ulations, so we turn to the thermal temperature of the 
gas for a better observational diagnostic. 

2.3.2. The M-T relation 

The temperature-mass relation for our models is shown 
in Figure ||. The temperature is a bolometric- luminosity 
weighted average across the cluster. In contrast to the M- 
o~ relation, even the most poorly resolved clusters follow 
the virial relation, with the exception of the open model, 
in which very low mass clusters have a higher temperature 
than predicted. This may be due to the earlier formation 
times for these clusters, as compared to the other models 



( [Navarro, Frenk fc White 1996Q . 

As before, the virial relations are fit with the 25 most 
massive clusters at z = 0, and the results, as parame- 
terized by fr are given in Table @ and plotted against 
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Fig. 4. — The virial mass and emissivity-weighted temperature of clusters for the same models and redshifts as in 
Figure ||. Lines are the virial scaling relations from equation (0), with normalizations given in Table ^. 
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cell size in Figure | There is no discernible trend with 
resolution. 

In the same Table, we also show the value of fx as de- 
termined by other simulations, all of which used smoothed 
particle hydrodynamics (SPH). The cluster simulations of 
Evrard (1991) produced a normalization very close to that 
obtained here, although he used a mass-weighted temper- 
ature. The six clusters of NFW again indicate a slightly 
higher value, although still within 10%. Recently, Evrard, 
Metzler & Navarro (1996, EMN96) examined simulations 
of 58 clusters (a sample which included the six of NFW) 
for a variety of models, including those with low Q. For 
A c = 250, they find a mean value of fx — 0.92. This 
is 17% above our mean result. A careful examination 
of their Figure 4 indicates that there are slight system- 
atic differences between the simulations, with clusters in 
flat universes implying somewhat lower values for fx , and 
SI = 0.2 (flat and open) models indicating higher nor- 
malizations. They also included a set of clusters which 
included substantial galactic winds, which appears to in- 
crease the virial normalization only very slightly. In it is 
interesting to note that the other physical contribution 
which we have not included here, cooling flows, would 
tend to decrease fx, although the expected effect is quite 
small. 

The range exhibited in Table ^ is surprisingly low, con- 
sidering the range of models, resolutions, numerical tech- 
niques and researchers who have examined this robust 
relation. The size of the scatter is also quite low, around 
10% (corresponding to a 15% scatter in mass). This im- 
plies that the temperature is an excellent virial mass indi- 



cator, even better than using /3-modcl estimates (Evrard, 
Me tier fc Navarro 199(f ) . 



2.3.3. The T-a relation 

The relation between temperature and velocity disper- 
sion is plotted in Figure [| along with equation (0) for 
three values of f3 (0.8,1,1.2). All simulations show curves 
which are steeper than predicted and there is a trend 
with mass, moving from low to high to values of [3. This 
is probably a reflection of the fact that low mass clusters 
are more poorly resolved and so have a lower a than pre- 
dicted (see also the M-a relation) . The best resolved clus- 
ters, in fact, have values above one. This follows from the 
M-T and M-a relations, since (3 = flj f T ~ 1-06. This 
value is close to the mean of NFW (1-07) but somewhat 
lower than the 1.18 found by Evrard (1991). Our result 
should most likely be viewed as a lower limit, since Fig- 
ure ^ indicates that the a values have not yet converged. 
The most likely explanation for /3 > 1 is incomplete ther- 



malization in the gas, although since the two value are 
computed with different weights, this is not completely 
straightforward. 

Observationally, the quantity f3 has been the subject of 
some investigation, with values being generally compati- 
ble with one (Girardi et al. 1996; Lubin & Bahcall 1992), 
although this is for the galaxy velocity dispersion rather 
than the collisionless component. See also Lubin et al. 
(1996). 

2.3.4- The M-Lboi relation 

Figure ^ shows the mass-bolometric luminosity relation 
for our three models. Shown are both the virial predic- 
tions from equation (|l0| ) as thin lines, and the resolution- 
adjusted equation (|15|), as thick curves. The latter rela- 
tion fits the simulated clusters quite well. As we demon- 



strated in section 2.2, a cluster's luminosity is diminished 



by a factor which depends on its mass (and the resolu- 
tion of the simulation). Another way to say this is that 
the core radius of the cluster is set by the the cell size, 
rather than scaling with cluster mass and density. Al- 
though the slopes differ, the variation in redshift for the 
two relations is very similar. Also, we should note that 
the curves are not fit (as for the M-T and M-a relations), 
but all have the same, somewhat arbitrary, normalization 
given in section [Q . 

2.3.5. The L{, i-T relation 

The final comparison is the bolometric luminosity- 
temperature relation (Figure @). Although it is not in- 
dependent of the previous four, it is the one with the 
best observational constraints. The observations indicate 
a slope of L ~ T 2 7 " 3 5 ( |Edge fc Stewart 1991a| ,b; [David 
t al. 1993). Although the slope of the simulated clus- 



ters is closer to the observations than that predicted by 
the equation ( |li"| ) , we argue that this is an artifact of our 
fixed numerical resolution. This result applies to Eule- 
rian simulations in general. Bryan et al. (1994b) found 
a slope in agreement with observations for the 2-10 keV 
luminosity-temperature (as opposed to the bolometric re- 
sult presented here). As we see here, this was partly 
due to resolution effects (but also resulted from the fi- 
nite band-pass which strongly curtails the luminosity for 
clusters under a few keV). Thus, the basic discrepancy 
between the observed and predicted bolometric L — T re- 
lation remains unresolved. Additional physics, such as 
galactic feedback, seems to be required (see also Evrard 
& Henry (1991) and NFW). 
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Fig. 5. — The velocity-dispersion and emissivity- weighted temperature of clusters identified for the same models and 
redshifts as in Figure 0, along with the scaling relation from equation (m). 
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Fig. 6. — The mass and bolometric luminosity of clusters identified for the same models and redshifts as in Figure ||. 
Thin curves are the virial scaling relation from equation (10), while thick lines arc from the resolution-adjusted version, 
equation (111). 
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Fig. 7. — The temperature and bolometric luminosity of clusters identified for the same models, redshifts and virial 
relations as in Figure B. The dot-dashed line indicates the observed relation L ~ T 3 at z = 0. 
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2.3.6. Scatter in the virial relations 

The scatter of clusters around the virial relations can 
be important because it provides a limit on how accu- 
rately masses can be determined, even in the absence of 
observational uncertainties. This is a factor in comput- 
ing distribution functions of observable quantities, such 
as temperature or luminosity, since a given mass is as- 
sociated with a range of temperatures. More concretely, 
given a mass distribution dn/dM, the resulting temper- 
ature function dn/dT (computed by multiplying by the 
virial relation dM/dT) must be convolved with the scat- 
ter. 

In Figure ||, we show cumulative distribution func- 
tions for our three canonical simulations. Although the 
fa- relation is somewhat tighter than fx, their accuracy 
in predicting the cluster mass is nearly identical (since 
M ~ er 3 ~ T 3 / 2 ). As we noted previously, the mean 
of the velocity dispersion distribution appears to shift 
slightly with resolution, but the temperature does not. 
The shape of the distributions are well described by a 
Gaussian (as verified by a Kolmogorov-Smirnov test) . 

In order to compare the size of the fluctuations to 
those observed in real clusters, we must use the L-T re- 
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Fig. 8. — Cumulative distributions of f a (top) and fx 
(bottom) for the 25 most massive clusters in our three 
primary simulations. The solid lines are Gaussian distri- 
butions fitted to the combined sample. 



lation. Although the slope of our predicted relation dis- 
agrees with that observed, we can compare the scatter 
around the two curves. Since real cluster temperatures 
also have observational uncertainties, we add the intrin- 
sic scatter and observed uncertainty in quadrature. The 
simulated clusters have a scatter of 11.4l 2 ' 3 %. The flux- 
limited sample of Henry & Arnaud (1991) gives 19.6^ 2 ' 7 % 
(similar to, but slightly higher than, the 16.4% found by 
Pen, 1996, for the same sample), while the larger sample 
of David et al. (1993) requires a smaller intrinsic scat- 
ter: 15.3±iJ%. These values are determined by fitting 
the L-T relation and adjusting <tt — fT until the result- 
ing x 2 indicates a good fit; uncertainties are one sigma. 
Although the scatter in the simulated clusters appears 
to be lower than that observed, there are two factors 
which would lead us to overestimate the observed value. 
The first, which is unlikely to be important, is ignoring 
the small observational uncertainties in the luminosity. 
More seriously, we have implicitly assumed that the tem- 
perature uncertainties are normally distributed. Since 
most realistic distributions have tail probabilities that are 
larger than Gaussian, this causes us to overvalue outliers, 
forcing the required intrinsic scatter to rise. In any case, 
the difference is not large, and so we tentatively conclude 
that the simulations reproduce the observed scatter. 

3. Distribution Functions 

The scaling relations discussed in the previous section 
can be combined with a prescription for computing the 
mass distribution function of virialized objects to make 
predictions of more easily observed distribution functions. 
Although there are many prescriptions for obtaining the 



linear initial perturbations (e.g. Lacey & Cole 1993; Ki- 


tayama & Suto 1996; 


Bond & Myers 1996a 


) , we adopt the 


Press-Schechter (PS) formalism ( 


Press & Schecter 1974) 



because it is fairly accurate, relatively simple, and widely 
used. 

The comoving number of virialized objects of mass M 
in mass interval dM is given by, 



dn 
dM 



2\ 1/2 p 



da 



M a{M) dM 



-vlll 



(18) 



where p is the mean density and v c — S c /a(M). The 
linear rms density fluctuations of the power spectrum is 
given by 

<t 2 (M) = [ W%(kR)P(k)d 3 k (19) 
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on the scale M — 4irR 3 p/3 with a top- hat smoothing fil- 
ter Wt(x) — 3(sin(x)/x — cos(a;))/x 2 [note the difference 
between the velocity dispersion a and the rms density 
fluctuation <r(M)]. The numerical value of 5 C has been 
the topic of some debate. The spherical collapse model 
gives S c as a slowing varying function of cosmology which 
is well fit by: S c = 1.686ft(z) c , where c = 0.055 (0.018) 
for flat (zero cosmological constant) models. However, 
given the shortcomings of the theory's theoretical under- 



pinnings (Bond et al. 1991), it is probably a better ap- 
proach to take S c as a free parameter and fit for it. 

In Figure ^ we show the result of comparing halos 
from the simulation (identified with the spherical overden- 
sity scheme described earlier) with predictions from equa- 
tion @ for the CDM270, CHDM512 and OCDM256 
models at a few different redshifts. The agreement is 
generally good, although there are systematically fewer 
low mass clusters in the simulation than predicted. This 
is due in part to the fact that those clusters are poorly 
resolved, but als o seems to be a gener al shortcoming of 
the PS method ( [Bond fc Myers 1996b| ). All the simula- 
tions produced fewer high-mass clusters with increasing 
redshift, although the evolution is strongest in the model 
with massive neutrinos and weakest in the open cosmol- 
ogy. The number of massive clusters at higher redshifts 
(particularly around z = 1) for the OCDM256 model is 
underpredicted. If this result is correct, then it makes the 
difference between open and flat models even larger than 
previously predicted by PS methods, however two com- 
ments are in order. First, this simulation has the smallest 
number of well-resolved clusters (CHDM512 has better 
resolution, while CDM270 has more and larger clusters). 
Second, the significance is small since the number of clus- 
ters involved is very small. Indeed, we should point out 
that these simulations are not the best to test the Press- 
Schecter methodology since purely N-body simulations 
can afford better resolution, and can use a larger box to 
produce more high mass clusters. Still, we press on to 
illustrate what we believe to be a more objective and in- 
formative way to constrain the parameter S c - 

We can be more quantitative about the fits by varying 
this parameter (fixed in Figure [9] by the spherical collapse 
prediction). Since the PS formalism is often used to con- 
strain theories of structure formation, knowledge of the 
uncertainty in its predictions is important. We construct 
a likelihood estimator by considering the distribution of 
clusters above a mass cutoff M* for each simulation. If 
we divide the range of masses into intervals AM that are 
sufficiently small so that the predicted number of objects 
in each range, n(M) — AMdn/dM <C 1, then the prob- 



ability of finding either or 1 objects in this range is 
given by the Poisson distribution: P(0) = exp[— n(M)] 
or P(l) = n(M)exp[-n(M)]. The relative likelihood of 
finding a particular set of clusters is the product of the 
probabilities for all the mass intervals, split into those 
with clusters and those without: P(5 C ) = f\ P(l) J] P (0). 
For AM <C 1, this can be written as: 



n 



dn(Mi 
dM 



■AM 



exp 



dn(M) 
dM 



dM 



(20) 

The product is over all the cluster masses M, above 
the cutoff M*; the distribution is normalized, so that 
J P(S c )dS c — 1. To imitate the simulation as closely as 
possible, we have set the power spectrum to zero for scales 
larger than the fundamental wavelength of the box. This 
has only a slight effect, shifting S c by about 1%. 

Figure [To] shows the resulting likelihood distributions 
for our three models. Since we are primarily interested in 
the higher mass clusters, we choose a mass cutoff (M*) of 
twice the non-linear mass (this is about 1 x 10 14 M Q for 
the OCDM256 and CHDM512 cases and 6 x 10 14 M Q for 
CDM270). The curves are reasonably well fit by Gaus- 
sians with means and one sigma uncertainties given by 
1.84 ± 0.12, 1.75 ± 0.12 and 1.72 ± 0.09 for CDM270, 
CHDM512 and OCDM256 respectively. These are all 
compatible with the spherical collapse prediction. 

The se values are in the range advocated by other au- 
thors (pfstathiou ct al. 1988fc [Narayan fc White 198S|; 



Carlberg & Couchman 1989: Bond et al. 1991), although 
some caution is in order since we are comparing results 
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Fig. 10. — The likelihood distributions for S c , as deter- 
mined by our three models at z — 0. 
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Fig. 9. — The mass distributions function for the same models and redshifts as in Figure g. Lines are Press-Schechter 
predictions with S c given by the spherical collapse model (see text). 



with different methods for identifying clusters and their 
masses. For more recent work which adopts the same 
methodology as employed here, the agreement is remark- 
able. Lacey & Cole (1996) used scale-free simulations 
with a variety of cluster analysis techniques. For the 
spherical ovcrdcnsity method with a top-hot filter, they 
find a best fit of S c = 1.96, however, this is over a con- 
siderable range in mass. For their rarest peaks (corre- 
sponding to the clusters discussed here), it appears that 
a normalization somewhat lower, but just above 1.68, is 
appropriate. Eke, Cole & Frenk (1996) found that a value 
around 1.75 gave good fits to very large (350 3 /i _1 Mpc 3 ) 
simulations for a variety of cosmologies. Borgani et al. 
(1997) argue for S c near 1.68 (for a top- hat filter) based 
on a series of simulations including a cosmological con- 
stant term. 

Finally, we note that a cluster identification based on 
mean overdensity does not necessarily produce a one-to- 
one correlation with bound halos. As we discuss in the 
section on temperature, this should not strongly affect our 
results, however it is important if one is interested in the 
details of bound halos. For a more complete discussion of 
the mass function of halos in the CDM model see Gelb & 
Bertschinger (1994). 

The high-end of the mass function is exponentially sen- 
sitive to the amplitude of the power spectrum, so it is 
worthwhile examining the non-linear mass {M n i), which 
controls the position of this cutoff. It is defined as the 
mass within a spherical volume for which the rms over- 



density from equation (|T^) is equal to S c . Figure |Tl| 
shows the evolution of this quantity, using both the lin- 
ear power spectrum and the fully non-linear (N-body) 
spectrum (for a more complete discussion in terms of the 
CDM model, see Jain & Bertschinger 1994). This figures 
shows both the increased normalization of the CDM270 
model, and the slower evolution of the OCDM256 mass 
function. Scale-free spectra with P(k) cx k n result in 
M n i cx (1 + z)- 6 /(™+ 3 ). The ordinate is log(l + z) so the 
these models would show up as straight lines on the plot. 
The linear results for CDM are closer to the non-linear 
than for the other models because of the shallower spec- 
tral profile (smaller logarithmic slope). In other words, 
the CHDM slope is close enough to P(k) ~ fc -3 that the 
extremely rapid evolution of M n \ requires more accuracy 
to obtain a good estimate of the non- linear mass. 

Although useful, the mass distribution is difficult to 
obtain from observations because of uncertainties in ob- 
taining the total mass out to the virial radius. Some- 
what easier is the temperature distribution. In Figure 
we plot this function for the canonical models as well as 
analytic predictions using Press-Schechter and the mass- 
temperature relation: equations (|l^) and (|J). From the 
results of section ||, we adopt the mass-temperature nor- 
malization fx = 0.77, producing reasonable agreement 
between the simulations and analytic results. Again, 
there is are strong difference between the evolutionary 
properties. When medium to high redshift temperature 
samples become available, they will provide strong con- 



16 



straints. 

We stress at this point that the temperature (and ve- 
locity dispersion) distribution functions are much more 
robust against variations in cluster identification shemcs 
than the mass function. This is true because the mass is 
a rising function of radius (since p ~ r~ 2 ), while, to first 
order, the temperature is constant. How can iV(> T) 
can be relatively independent of the cluster identification 
method, while N(> M) — from which it is derived — is 
not? In fact, this is only true if the N(> M) and the T-M 
relation are calibrated consistently. In other words, if we 
change the way we find the cluster mass (within reason), 
the T-M relation and the mass function will both change 
such that the resulting iV(> T) is unchanged. This fea- 
ture of the temperature distribution is particularly impor- 
tant when it comes time to compare with observations, 
which are rarely determined to the virial radius. 

Since luminosities are more easily obtained observa- 
tionally, we use equation ( [[o| ) to determine the bolomet- 
ric luminosity distribution with Bremsstrahlung radiation 
and zero metallicity (this constraint will be relaxed in the 
next section) . These are shown as thin lines in Figure |l3| . 
While the simulated clusters do not match these results, 
including the effects of finite resolution with equation ( |To| ) 
does produce agreement, at least for the relative bright 
clusters. The 'fully' resolved Press-Schechter predictions 



1(T 
10" 

10 13 
10 12 

*10 11 
f 10 1 ° 

10 9 

10 8 
10 7 




• • CDM270 (non-linear) 

c o CDM270 (linear) 

■ ■ CHDM512 (non-linear) 
0---0 CHDM512 (linear) 

♦ » OCDM270 (non-linear) 

-o OCDM270 (linear) 



0.1 



0.3 0.5 



1.0 



2.0 



Fig. 11. — The evolution of the non-linear mass (see text 
for definition) for our three primary models, computed 
both through the linear and non-linear (N-body) power 
spectra. 



show that, for the brightest clusters, there is positive evo- 
lution with redshift in the OCDM model but negative evo- 
lution for CHDM. However, the exact behaviour depends 
on the adopted L-T relation. Since the fixed resolution 
produces an L — T relation closer to that observed (albeit 
for numerical reasons), it may be the better predictor of 
evolution. In that case, the OCDM model is nearly con- 
stant over redshift, while the CHDM evolution is strong 
negative. 

Perhaps more importantly, the difference between the 
thin and bold lines indicates the kind of error engendered 
by our limited spatial resolution. Roughly speaking, there 
are two effects. The first is the obvious reduction in the 
luminosity of all clusters, shifting the curves to the left. 
The other is to add an additional component of negative 
redshift evolution, exaggerating the negative evolution at 
high luminosities and decreasing the positive evolution 
at the low end, causing the curves to tilt slightly. This 
second effect is relatively small, about the size of Poisson 
fluctuations due to the limited number of clusters. 

3.1. Metallicity Effects 

In this section we address more realistic X-ray emis- 
sivity effects: limited band-passes and line emission from 
multiply ionized heavy elements in the plasma: metallic- 
ity. As we will show, for moderate and high temperature 
clusters, the ratio of line to continuum emission is actu- 
ally quite small, only a few percent, making the accurate 
determination of a cluster's metallicity difficult. On the 
positive side, this make our results largely insensitive to 
the assumed metallicity. Nevertheless, the high resolution 
spectroscopy performed by the ASCA satellite, as well as 
work done with ROSAT and other missions, seems to in- 
dicate that a value of Z /Zp, = 0.3 in solar units is reason- 



;t al. 1994| ) 



able for ma ny clusters ( Kowalski et al. 1993 ; Markcvitch 



We use a Raymond-Smith code (1977, 1992 version) 
to compute emissivities which are shown in the top panel 
of Figure [li] for a variety of metallicities and band-passes 
as a function of temperature assuming a number density 
of electrons of 1CP 3 cm~ 3 . The effect of line emission 
from metals is small for T > 2 kcV. This means that 



previous results ( Bryan ct al. 1994a ,b) computed with the 
assumption of zero-metallicity will not be substantially 
affected in the region where observational results are most 
plentiful. Nevertheless, this will prove important for small 
clusters and groups of galaxies so we extend our previous 
analytic results to include this effect (although, of course, 
our simulations do not include many other effects that 
may prove important for small clusters such as stellar 
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Fig. 12. — The temperature distribution functions for the same models and redshifts as in Figure g. Lines are 
predictions from the extended PS model with S c given by the spherical collapse model and fx = 0.77. 




Fig. 13. — The bolometric luminosity distribution functions for the same models and redshifts as in Figure 0. The 
Press-Schechter results for two mass-luminosity relations are also shown: thin lines use equation ( |l0| ) while bold lines 
include the effects of limited resolution: equation (ffa). 
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feedback and radiative cooling). 

To do this, we fit simple expressions to the full result 
by ignoring the temperature variation of the gaunt factor 
and approximating the effect of metallicity as a power law 
with a cutoff. The resulting formula is given in relation 
to the bolometric luminosity {L meta i — e(Z, T, v\, v^L^oi-, 
where v\ and v-i specify the band-pass) by: 



e(Z, T, vi,u 2 ) 



exp(- 



(kT/2 keV) 
1 



kT ' 



exp( 



kT 

kT < 2 keV 
kT > 2 keV 



(21) 



The effect of metallicity depends on the bandpass used 
and cannot be simply modelled, so we parameterize it 
through 7 and tabulate values appropriate for a few com- 
mon choices: 7 = 0.6 for 0.1-2.4 keV; 0.9 for 0.5-2.4 keV; 
0.0 for 2-10 keV; and 0.7 for bolometric luminosities (we 
use gff = 1-3, and 24% He by mass). This fit is accurate 
to only about a factor of two, but will be sufficient for our 
purposes; a few examples are shown in the bottom panel 



of Figure 14 



With this extension we can now compute more ob- 
servationally oriented luminosity functions and compare 
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Fig. 14. — The X-ray emissivity as a function of temper- 
ature for two different metallicities and three different 
band-passes from a Raymond-Smith code (top panel). 
The relative X-ray emissivity (normalized to the bolo- 
metric emissivity at that temperature for a metal-free 
plasma) as thin lines (bottom panel) . Also shown in thick 
lines are the approximate fits from equation (21). 



them to the luminosity function computed directly from 
the simulation. The numerical cluster luminosities are 
computed on a cell-by-cell basis so they include the spa- 
tial variation of temperatures (see section |]), which the 
analytic models are missing. They also use emissivities 
calculated from the full Raymond-Smith code. Figure |l5| 
shows the luminosity function in the 0.1-2.4 keV band- 
pass for a representative case (CDM270 at z = 0). The 
PS model overpredicts the number of low- luminosity clus- 
ters, just as it does for low mass clusters, however the gen- 
eral effect of metallicity is included correctly. This also 
demonstrates that the effects of metallicity are negligible 
above about 10 44 erg/s and relatively slight above 10 43 
erg/s. 

4. Cluster structure 

In deriving some of the analytic results we assumed a 
spherically symmetric isothermal profile for both the gas 
and dark matter. In this section we examine the simu- 
lated clusters in order to determine the accuracy of these 
assumptions. A more complete analysis will be presented 
in Bryan & Norman (1997b). 

In Figure [l6] we show profiles of temperature and 
the one-dimensional velocity dispersion for the five most 
massive clusters in our two canonical models. These 
are normalized by their appropriate virial values (with 
fx = fa = !)■ To compute the profile, we redetermine 
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Fig. 15. — The 0.1-2.4 keV luminosity distribution func- 
tion with (dashed) and without (solid) metallicity for the 
CDM270 model. The extended Press-Schechter results 
are also shown, with the usual meaning for thin and thick 
lines. 
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the cluster centers by adopting the center of the cell with 
the highest gas density within r V i r /2 of the original center 
found through the iterated spherical overdensity method. 
This procedure does a good job of finding the core of the 
largest mass clump. The innermost point plotted is two 
cell widths from the cluster center. This is approximately 
our resolution limit. 

The spherically-averaged temperature profiles are com- 
patible between models and show a very slowly falling 
profile to about 0.7r V i r , a somewhat steeper slope (r -0 - 5 ) 
to about two times the virial radius and then a very sharp 
fall off beyond that. The velocity profiles are close to 
flat, although they appear to show a slight dip at and 
just before the virial radius. They are also slightly below 
their respective a V i r values (implying f a < 1). However, 
this tendency may be declining as we resolve further to- 
wards the center. The hot component (we plot only the 
cold particles in Figure |l6|) follows almost exactly the 
same profile at large r, but is systematically higher at low 
radii. The temperature profile shows much greater vari- 
ation and is less compatible with an isothermal model, 
even after spherical averaging. The temperature and a 
profiles are in agreement with the factors fx ~ 0.8 and 
f (j ~ 0.85-0.9 derived earlier. 
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Fig. 16. — Temperature (top panel) and dark matter one- 
dimensional velocity dispersion (bottom panel) profiles 
for the CDM270 and the CHDM512 models at z = 0, 
normalized by their respective virial values. 



5. Conclusion 

In this paper, we have shown that X-ray clusters pro- 
duced in Eulerian simulations agree well with the scaling 
relations involving mass, temperature and the collision- 
less velocity dispersion. The luminosity behaviour is more 
affected by resolution but can be simply and accurately 
modelled. The predicted bolometric L-T relation does not 
match that observed. We argue that this is mostly the 
fault of the luminosity prediction since it is very sensitive 
to the structure of the clusters, while the temperature is 
not. Indeed, the mass-temperature relation appears to 
be remarkably robust and, as we have demonstrated, de- 
pends very little on the input physics, numerical method, 
resolution, or cosmology. Our result, combined with a 
review of the available literature, indicates that fx ~ 0.8 
in the notation of equation (|^). We also demonstrated 
that the isothermal profile assumed in computing the nor- 
malization of the scaling relation between M and a is a 
reasonable approximation for the dark matter (although 
there is evidence for f a < 1). 

We stress here the difference between fx, fa and /3. 
The first two adjust the virial scaling relation normaliza- 
tions as compared to the hydrostatic isothermal sphere 
assumption, while f3 is the ratio of the collisionless 'tem- 
perature' to the gas temperature, equation (]?]). These 
three quantities are related through (3 = fH fx- Since 
prescriptions such as the Press-Schechter formalism pre- 
dict the distribution of masses, the important normal- 
izations are f a and fx (of the M-a and M-T relations, 
respectively), and not f3 , as is occasionally assumed (e.g. 
Eke, Cole k Frenk 1996| ). 



Based on these findings and the Press-Schechter pre- 
scription for computing the differential number density 
of virialized halos, we computed the temperature dis- 
tributions functions which are in good agreement with 
those derived from the numerical simulations (over their 
range of validity). In Figure [l7|, we show the tempera- 
ture functions for our three models along with the 95% 
uncertainty in 5 C . The observations (Henry k Arnaud 
199l|; |Eke, Cole k Frenk 1996| ) are also shown, com- 



puted as N(> T) = J2xi>T l/ v max,i (where V m ax,i is 
the maximum volume for which the cluster with temper- 
ature Ti and flux fi could have been observed). This 
shows that the CDM270 model predicts substantially too 
many clusters, while the CHDM512 and OCDM256 mod- 
els are (very) marginally in agreement. The fits would 
be much improved if the power spectra were renormal- 
ized to cr 8 = 0.55,0.63,0.85 for CDM270, CHDM512 & 
OCDM256. To retain agreement with COBE, the CHDM 
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model would require only a slight tilt. This result dif- 
fers slightly from Ma (1996) due to the addition of a 
second neutrino. Otherwise, for the same Q rm s~ps we 



would have to increase a a by 10% (Stompor, Gorski & 
Bafrday lgglh. The OCDM270 model could also be ad- 



justed slightly to match both COBE and clusters, how- 
ever, CDM would require a large tilt in the spectrum 
(n ~ 0.8). 

Analytic luminosity functions computed with the ef- 
fects of finite resolution (Figure |l|) agree well with the 
simulations and allow us to gauge the impact of resolu- 
tion as a function of luminosity and redshift. An exten- 
sion which accounts for metallicity and limited bandpass, 
equation (pl|), can be used to compare directly against 
observational data. Just as we used the temperature 
distribution function, we could also use the luminosity 
to compare to observations, however, the observational 
luminosity-temperature relation does not agree with the 
scaling relations. Using the observed L-T relation to com- 
pute n(> L) from n(> T) does not add any additional in- 
formation, since if we match the temperature function we 
must also match the luminosity function. However, recent 
observations with the ASCA satellite have produced a 



sample of temperatures at moderate redshift ( Mushotzky 



& ' charf 1997). The sample is not complete so we cannot 
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Fig. 17. — The PS temperature distribution functions 
for our three primary models at z — using the best 
fit S c for each model. Shaded regions are the two-sigma 
uncertainties in S c . Observations from Henry & Arnaud 
(1991) are shown as a solid line. 



use it to construct a temperature function, however it can 
be used to constrain evolution in the observed L-T rela- 
tion. Combined with the method for predicting n(> T) 
described in this paper and the upcoming high redshift lu- 



minosity samples (Romer ct al. 1997; Scharf et al. 1997), 
we will be able to place much tighter constraints on the 
models based on cluster evolution. 

This work is done under the auspices of the Grand 
Challenge Cosmology Consortium and supported in part 
by NSF grants ASC-9318185 and NASA Long Term As- 
trophysics grant NAGW-3152. 
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